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Continuous-Time Quantum Monte Carlo Approach to Singlet-Triplet Kondo Systems 

Shintaro Hoshino*, Junya Otsuki and Yoshio Kuramoto 

Department of Physics, Tohoku University, Sendai 980-8578 
(Received April 10, 2009) 

Dynamical properties are studied numerically for a variant of the Kondo model with singlet 
and triplet crystalline electric field (CEF) levels where Kondo and CEF singlets compete for the 
ground state. Using the continuous-time quantum Monte Carlo method, we derive the t-matrix 
of conduction electrons and dynamical susceptibilities of local electrons without encountering 
the negative sign problem. When the CEF splitting is comparable to the Kondo temperature, 
the dynamical response has only a quasi-elastic peak. Nevertheless, the local single-particle 
spectrum shows an energy gap in strong contrast with the ordinary Kondo model. 
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1. Introduction 

In Pr and U compounds with even number of lo- 
calized / electrons per site, the crystalline electric field 
(CEF) ground state can be a singlet.^' Interestingly, 
some of these systems show the Kondo effect, while oth- 
ers do not. Hence it has long been a focus of interest 
how the competition works between the CEF singlet and 
the Kondo singlet under the exchange coupling with con- 
duction electrons. This model is called the singlet-triplet 
(ST) Kondo model in the rest of the paper. Theoretical 
methods for studying the ST Kondo model include the 
poor man's scaling^^ and the numerical renormalization 
group 

(NRG). 4-6) For Pr skutterudites, the non-crossmg 
approximation (NCA) and the NRG has recently been 
used for deriving the dynamics of the system.^' In or- 
der to clarify the details of competition not only for the 
ground state but also for finite temperatures, however, 
more accurate theory is desirable. 

In this paper, we derive the dynamical property of 
the ST Kondo model at finite temperature using a new 
accurate method. The competition between CEF and 
Kondo effects is studied by the continuous-time quan- 
tum Monte Carlo method (CT-QMC), which has been 
first proposed by Rubstov et al.^^ The original CT- 
QMC is formulated by perturbation expansion in terms 
of Coulomb interaction. Werner and Millis developed a 
way how the expansion is performed in terms of the 
hybridization between the impurity and the conduction 
electrons. ^''^ Recently, the CT-QMC method is extended 
to the Coqblin-Schrieffer model and the Kondo model 
by Otsuki et al.^^^ In this paper, we apply the CT- 
QMC to the ST Kondo model proposed in ref. 12 where 
the singlet-triplet states interact with the conduction- 
electron spins. It is found that the calculation can be 
performed without negative sign problem. Hence, the nu- 
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merical results shown in this paper are highly accurate. 
The single-particle excitations and the dynamical suscep- 
tibilities are obtained for the first time in the ST Kondo 
model. 

In §2, we introduce the ST Kondo model and classify 
the fixed points according to the signs of two exchange 
interactions. Then in §3 we formulate the CT-QMC for 
the ST Kondo model, and discuss the details of actual 
Monte Carlo simulations. We introduce correlation func- 
tions which give static and dynamic susceptibilities in 
§4. Numerical results for the ST Kondo model are pre- 
sented in §5 for the doublet ground state, and in §6 for 
the singlet ground state. We summarize our results in §7. 
In Appendix, we discuss technical aspects of extending 
the CT-QMC to the ferromagnetic exchange. 

2. Singlet-Triplet Kondo Model and Its Fixed 
Points 

We consider the CEF singlet-triplet system interact- 
ing with conduction electrons. The singlet-triplet levels 
can be represented as a spin dimer with two pseudo spins 
S'l and S2 - In the pseudo-spin representation of the CEF 
levels, the spin singlet and triplet describe the CEF sin- 
glet and triplet, respectively. The CEF splitting is rep- 
resented by AcEF- The ST Kondo model is then written 



as 



12) 



H 



her 



2{JiSi + J2S2) ■ Sc 



iCEF 



Si ■ S2 



(1) 



where ^fc is the energy of the conduction electrons mea- 
sured from the chemical potential, and the operator 

1 



(2) 



describes the conduction-electron spin at the origin. Here 
Cct(cJ.) is an annihilation (creation) operator of the con- 
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duction electron. 

Using operators X* ~ Si + S2 and X" = Si — 
S2, the interaction term between the pseudo spins and 
conduction spin in eq.(l) is rewritten as 

Wint = (/tX'+/,X^).Se, (3) 

where It = Ji + J2, Is — Ji — J2- The operator X* de- 
scribes transition within triplet states, and connects 
a singlet state with triplet states. In this paper, we 
use the rectangular density of states for the conduction 
electrons: 

p{e) - pce{D - (4) 

where I? is a band width and pc ~ 1/2D. We use the 
unit Z) = 1 in numerical calculations. 

In the CT-QMC method, it is difficult to treat the 
CEF term. Hence we impose the condition Aqef = 0. 
Even under this condition, an effective CEF splitting 
arises since two pseudo spins interact with each other 
through the conduction electrons. Let us derive the ef- 
fective interaction between the pseudo spins by applying 
the second-order perturbation theory to the Hamiltonian 
(1) with AcEF = 0. Taking the expectation value with 
respect to conduction electrons with the density of state 
(4), the effective CEF splitting in the T ^ limit is 
given by 

In 4 

AcEF = -^yJiJ2- (5) 

If Ji and J2 have different signs, Acef is positive and 
the effective interaction stabilizes the CEF singlet. 

We classify the parameter space into (I) Ji , J2 > 
0, (II) Ji > 0, J2 < and (III) Ji, J2 < 0. Figure 1 
shows a schematic phase diagram for Acef = 0. The 
ground state of each region is (I) doublet, (II) singlet and 
(III) triplet. The residual entropy in (I) is understood 
as due to a remaining free spin after the Kondo effect 
compensates one of the two spins. On the other hand, 
the entropies in (II) and (III) correspond to CEF singlet 
and triplet, respectively. 

If a coupling constant increases toward a posi- 
tive direction, the Kondo effect for the spin S^ becomes 
stronger. If both | Ji| and | J2I increase, the effective CEF 
splitting also increases according to eq.(5). Hence we can 
see the competition between Kondo and CEF effects by 
studying representative points of (Ji, J2) in Fig.l. Ac- 
cording to NRG,^'^-' the boundary between singlet and 
doublet is in fact in the region (II) close to the Ji axis. 

In the region (II), the character of the singlet changes 
from CEF to Kondo as ( Ji, J2) goes away from the ori- 
gin. From the scaling theory for the Kondo model, we 
estimate the Kondo temperature as 

which contains only the antifcrromagnetic coupling Ji. 
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Fig. 1. Fixed points of the ST Kondo model for Acef = 0- 
The properties of this system are symmetric with respect to 
the straight line Ji = J2- The parameters indicated by circles, 
squares and triangles are used later in numerical calculation. 



Since the ferromagnetic interaction is renormalized to 
in the Kondo model, we consider eq.(6) as the energy 
scale of the Kondo effect. On the other hand, the char- 
acteristic energy for the CEF effect is given by eq. (5). 
Using these two characteristic energy values, we define 
the ratio: 

r^J^. (7) 
Acef 

The Kondo effect is negligible for r ^ 1 , and competition 
against the CEF effect arises for r ^ 1. 

3. Monte Carlo Procedure 

The algorithm for the Kondo model has been devel- 
oped in ref.ll. Here we describe how to deal with the ST 
Kondo model with Acef = using CT-QMC. Since the 
original algorithm is designed for the antiferromagnetic 
exchange, we have extended the algorithm so as to be 
applicable also to the ferromagnetic exchange as shown 
in Appendix. 

We rewrite the Hamiltonian using permutation op- 
erators as follows: 

H = Hc+ Hint , (8) 

= XI ^kci^Cka - '^^-J^ clca, (9) 
ka (7 

Hint = X! •^A'X^^^'T'lCa'CtT -aM^(T<T') 

/i— 1.2 aa' fJ- 

(10) 

which form is suitable to apply the CT-QMC. Here X^^, 
is the operator that changes the p-th pseudo-spin state 
from a' to a. The parameter is introduced to avoid 
the negative sign configuration (see Appendix for detail) . 
We take to be for the ferromagnetic coupling and 
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1 for the antiferromagnetic coupling. The constant term 
J2fj, nis-y be neglected. The partition function Z = 

Tj.g-/3-H jg factorized as Z = ZiZ, with Z^ = Tr^ e"^^-. 
We obtain Zi as 



Vq W{q), 



ordered 

dn-. 



W{q) = 



it^k} {cTk} Wk} 



(11) 



(12) 



nTr[T.<,,(ri)---<^,, (r,J], (13) 



(fc,K},{Mfc},K},K.}), 



(14) 



where 'ordered' in eq.(12) means that the configuration 
{Tfc} is aligned as /3 > > • ■ • > r2 > ti > 0. The nota- 
tion I means the interaction picture for the time evolu- 
tion. The order k of the expansion consists of components 
ka- and fc^ with the relation k^ = fc^ = fc. The 
suflSces (T and /i correspond to a set of operators cJ.C(j 
and to X'', respectively. 

The conduction-electron operators are grouped by 
their spin indices. Using Wick's theorem, the spin a 
conduction-electron part is represented as determinant 

(k ) 

of ka- X ka matrix " whose elements are composed by 
the free conduction-electron Green function g{T). A ran- 
dom walk in the configuration space {g} with the weight 
W{q) enables us to perform the integral (12). When k^ 
increases, it takes much time to update the matrices of 
the conduction-electron parts in numerical calculations. 
Note that the most important value of ka depends on the 
strength of the interaction J^. Since two pseudo spins in- 
teract with conduction electrons, the calculation becomes 
heavier than the Kondo model, which has only one local 
spin. We also note that the antiferromagnetic coupling 
makes ka larger than the ferromagnetic coupling. 

We represent the configuration q by a diagram as in 
ref. 11. Figure 2 compares the diagram for the Kondo 
model and the ST Kondo model. The ST Kondo model 
has two pseudo-spin components labeled by 1 and 2. In 
order to change a configuration in Fig. 2(a), we add a 
new operator Xaa' to the local spin configuration, and 
c^^,Ca to the conduction parts. In Fig. 2(b), on the other 
hand, we first choose either /i = 1 or 2 and then add 
operators as in Fig. 2(a). Thus, the updating process in 
the ST Kondo model is the same as in the Kondo model 
except for additional choice of a pseudo-spin component. 

In the simulation of the ST Kondo model, we have 
observed negative weight configurations at the rate of 
about 10^^ at low temperature, which is to be compared 
with 10~^ in the Kondo model. The occurrence here is 
largely due to rounding errors, since it depends on the 




Fig. 2. Diagram of a configuration q with k\ = 2, k2 = 1, fc| = 2, 
fc| = 1 for (a) the Kondo model and (b) the ST Kondo model. 
The X-operators act on the local state, while the incoming and 
outgoing arrows denote the annihilation and creation of conduc- 
tion electrons, respectively. 



numerical treatment of the Green-function determinant 
(k ) 

det D'a "'. The weight is small enough to perform the 
simulation accurately. 

Figure 3 shows the probability distribution with re- 
spect to k^ and ka defined by eq.(13). In the region (I), 
as shown in (a) with (Ji, J2) = (0.3, 0.3), the distribution 
of ka has two peaks, which are caused by spin fluctua- 
tion; if an |-spin is on the left peak in a snapshot of the 
Monte Carlo simulation, another spin with | is on the 
right peak. At low temperature, these two peaks tend to 
separate completely. In this case, we cannot move around 
all the configuration space {q}, since each spin is trapped 
on the different peak during the simulation. We can re- 
solve this problem by introducing a sweeping procedure 
that flips all the pseudo-spin states. This update can be 
done safely because it does not change the value of W{q) 
in the ST Kondo model. On the other hand in the re- 
gion (II), as shown in (b) with (Ji, J2) = (0.3, —0.3), the 
probability distribution has only a single peak. Hence we 
do not need the sweeping procedure. 

In the region (III), as shown in (c) with (Ji, J2) = 
(—0.3, —0.3), the probability distribution with respect to 
ka has triple peaks due to fluctuation among three com- 
ponents. The number of degrees of freedom means that 
the triplet states are stabilized against the singlet state 
by the effective CEF splitting. In the Monte Carlo sim- 
ulation, presence of an f-spin on the left peak means 
presence of a |-spin on the right peak, which is the same 
as in the region (I). On the other hand, if an |-spin is 
on the central peak, then another spin with | is also on 
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Fig. 3. The probability distribution P(fc) at T = 0.01 in (a) re- 
gion (I), (b) region (II), and (c) region (III). The exchange pa- 
rameters are shown as crosses in Fig. 1. 



the central peak. At low temperature, these three peaks 
separate completely. We cannot move around all the con- 
figuration space even by the sweeping spin flip because 
of the central peak. For a simulation in the region (III) , 
therefore, wc need another procedure to mix the central 
peak and the side peaks. Physically, the region (III) is 
less interesting because both exchange couplings Ji , J2 
renormalize to zero. Hence we do not consider this re- 
gion in the rest of the paper. 

In the Monte Carlo simulation, we estimate statis- 
tical errors from 20 bins of data. To obtain dynamical 
quantities on the real-frequency axis, we have performed 



analytic continuation using the Pade approximation. Al- 
though the Pade approximation does not take statistical 
errors into account, the data of the CT-QMC are accu- 
rate enough to obtain reliable dynamics. This aspect has 
already been demonstrated in ref.ll. In particular, at low 
temperature the Pade approximation well reproduces the 
spectrum, since the data interval 27rT on the imaginary 
axis becomes narrow. 

The ST Kondo model has the particle-hole symme- 
try in the single-particle spectrum. It then follows that 
the t-matrix in the Matsubara frequency domain should 
be pure imaginary. We have imposed the condition be- 
fore analytic continuation, by neglecting tiny real part 
arising from statistical errors. As a result, the spectrum 
maintains the particle-hole symmetry. 

4. Correlation Functions 

As the most fundamental quantity in the system, 
the imaginary time correlation function is discussed now. 
The correlation function has a label of pseudo spins, and 
is defined by 

(15) 



where n^cr = Xj^^^ is the number operator of the 11- 
th pseudo spin with a spin ct, and tilde means devia- 
tion from the mean number: n^o- = "-/^cr — {n^a)- Sim- 
ilarly we define a pseudo-spin correlation function by 
xZi'^) = {SI,{t)SI) where S*^ = Y.^^n^„ is the /x-th 
pseudo-spin magnetic moment [a — ±1/2). We intro- 
duce a correlation function Xm in terms of x^a' ^ 



C 



Acrcr 



ir)], 



(16) 



where C = 1/4 is the Curie constant, and a = —a. The 
correlation function has the symmetry Xm ^ x'm- The 
static susceptibility is obtained by integrating the corre- 
lation function from to (3. With use of X' and X^, we 
define related correlation functions as follows: 



Xt,s(T) = {Xr{T)Xr) = Xuir) + Xuir) ± 2xJ.?(r). 



(17) 



The physical magnetic moment is given by = 

a^S^ where the coefficient depends on wave func- 
tion of the local states. Then the magnetic correlation 
function xj is represented by xj(t) = J^ni, 
In this paper, we show results for susceptibilities Xm 
separately in order to see responses of each pseudo spin. 
Other multipoles like a quadrupole can be also written 
in terms of Si and 82- 

The projection operator Pg onto the singlet state is 
given by Pg = —Si ■ S2 + 1/4, and the triplet projection 
is given by Pt = 1 — -Ps- We can derive the singlet occu- 
pation rate with use of the correlation function In 
the isotropic system, the singlet occupation rate is given 
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Fig. 4. Singlet occupation rate (Ps) as a function of temperature. 
The parameters (Ji, J2) are shown as circles and triangles in Fig. 
1. 

by 

(^^s> = i-3x3i(r = 0). (18) 

In the high-temperature hmit, we obtain (Pg) = 1/4 from 
eq.(18), since there is no correlation between different 
pseudo spins in this limit. 

It is convenient to characterize the CEF states in 
terms of the singlet occupation rate {Ps)- Figure 4 shows 
(Pg) for the region (I) and (II), with contrasting depen- 
dence on temperature. The value of (P,) changes most 
significantly at temperatures corresponding to the the ef- 
fective CEF splitting. In the case of ( Ji, J2) = (0.2, 0.2), 
the singlet occupation rate tends to due to the rela- 
tion /g = Ji — J2 = 0. Namely, only the triplet gains the 
interaction energy which is always negative in second- 
order. On the other hand, with (Ji, J2) = (0.2, 0.05), the 
singlet state also participate in the ground state to give 

(Ps)7^0. 

In the region (II), (Pg) comes close to unity with 
r -C 1 as defined by eq.(7). For example, in the case of 
(Ji, J2) = (0.1,-0.1) we obtain r = 0.0033 « 1, and 
the CEF effect is dominant over the Kondo effect. For 
larger Ji = — J2 in the region (II), the Kondo effect be- 
comes important, and the singlet occupation rate (Ps) 
decreases. For example, the case (Ji, J2) = (0.3,-0.3) 
gives r = 0.29 ~ 1, and (Pg) does not tend to unity even 
at low temperature. 

In the following, we discuss the regions (I) and (II) 
in more detail. 

5. Doublet Ground State by Underscreened 
Kondo Effect 

Let us first concentrate on the region (I) where the 

underscreened Kondo effect occurs. As a typical example, 
we take the case Ji = J2 = J > Q where the interaction 



Hamiltonian is written as 

Wint - J{Sl + S2) ■ Sc. (19) 

In the strong-coupling limit, we can neglect the kinetic 
energy of conduction electrons. Then the ground state 
\ga) is a doublet composed by linear combination of |5f + 
5*1, Sg) where only the triplet {S = 1) part enters in 
Sf + 82- Namely we obtain 




and the time-reversal partner |g ].) with the correspond- 
ing energy Eg = —J. We have the relation Xm — Xm 
because of the condition Ji = J2. The squared effective 
moment is given by 

TXt = |(g T \St + Sl\g t)|2 = 4/9 ~ 0.44 (21) 

for the ground state. 

Figure 5 shows the numerical results for 
the region (I). The property Xm > ^ indicates the fer- 
romagnetic correlation between the pseudo spins Si and 
82- At high temperature, the result tends to the Curie 
law Txm /C = 1- The effective moments seem to become 
constant at sufficiently low temperature. Note that the 
value is unity if we have a free spin with 5 = 1/2. In Fig. 
5(a) with Ji = J2 = 0.2, both Txu and Txu tend to the 
same value in the low temperature limit. This is because 
the condition Xm = Xm = Xt/4 is satisfied at sufficiently 
low temperature due to the absense of Ig = Ji — J2. If 
there is no Kondo effect, we expect Txt = 2/3 ~ 0.67 as- 
sociated with the triplet ground state. On the other hand, 
the strong coupling limit gives Txt ~ 0.44 as given by 
eq.(21). The computed value Txt — 0.4 in the low tem- 
perature limit is close to the strong coupling limit, and is 
far from the values for a free spin with either S ~ 1/2 or 
S = 1. On the other hand, if we have Ji ^ J2 as shown 
in Fig. 5(b), we no longer have the condition Xm ~ Xm- 
In this case, the effective moment tends to another value 
which depends on the ratio Ji/J2- 

The single-particle spcctriim of the underscreened 
Kondo system is also interesting but has not been inves- 
tigated so far to our knowledge. The relevant quantity 
is — Imfo-('^ + ^S) where ta-{z) is the impurity i-matrix 
of spin a. We derive —lmta{uj + iS) following the proce- 
dure of ref.ll. Figure 6 shows the result at low temper- 
ature. We plot two lines corresponding to if and for 
each parameter set. The coincidence of the two indicates 
that the Fade approximation is reliable. There appears 
a peak at the Fermi level, which is due to the under- 
screened Kondo effect. Note that the spectral shape is 
far from Lorentzian, but is characterized by two differ- 
ent energy scales, which is most evident in the case of 
{•-h- J2) = (0.2,0.05). This feature is in marked contrast 
with the spectrum in the ordinary Kondo model, which 
has a single energy scale Tk- 



5/9 



J. Phys. Soc. Jpn. 



Full Paper 



s 





0.001 



Fig. 5. Tx'^ /C as a function of temperature in the region (I). 
The symbols 11, 22 and 12 mean Xmi Xm ^'^'^ respectively. 
Parameters are shown as circles in Fig. 1. 



6. CEF Singlet vs Kondo Singlet 

We now discuss physical properties in the region (II) 
where the ground state is a singlet state. The competition 
between the Kondo effect and the CEF effect determines 
the behavior of this system. 

6.1 Impurity t-matrix 

In the perturbation theory with respect to Ji and J2 , 
the imaginary part of the i-matrix has a threshold singu- 
larity since the conduction electron whose energy is lower 
than the effective CEF splitting cannot be scattered by 
the local state. 

Figure 7 shows numerical results for —lmt„{LO + i5) 
in the region (II), obtained by the procedure of ref.ll. 
Figure 7(a) shows the case with a gap structure. In the 
case of a small value of r defined by eq.(7), the spectral 
function is similar to the step function. This is because 
the CEF effect is dominant compared to the Kondo ef- 
fect. A typical case with (Ji, J2) = (0.1, —0.1) (r <C 1) is 
shown also in the inset of Fig. 7(a). As r grows, the gap 
structure becomes obscure. 

We now fix Ji = 0.3 and vary J2 from zero to — Ji. 
Figure 7(b) shows the resuhs. With ( Ji, J2) = (0.3, -0.3) 
(r ~ 1), for example, the spectral shape has the two 



E 
T 




0.05 0.1 



Fig. 6. t-matrix at T = 0.001 in the region (I). Two lines corre- 
sponding to up and down spins are given for each parameter set 
(Jl, Ji). The good agreement of the two lines indicates the accu- 
racy of Fade approximation. The parameters (Ji, J2) are shown 
as circle points in Fig. 1. 



broad peaks at uj ^ ±0.1. As IJ2I decreases with Ji 
fixed, the effective CEF splitting decreases according to 
eq.(5) while the Kondo temperature Tk defined in eq.(6) 
is constant. Hence the ratio r increases with decreasing 
IJ2I. The gap in the spectrum becomes narrower as IJ2I 
decreases, reflecting the smaller effective CEF splitting. 
Note that the temperature T = 0.001 is not low enough 
compared to the effective CEF splitting for J2 = —0.05 
and —0.1. Hence the spectrum at the Fermi level is not 
zero. In the case of J2 = 0, which gives r — 00, these 
two peaks combine into one. The case corresponds to the 
ordinary Kondo effect, however with an extra free spin. 

6.2 Dynamical Susceptibility 

We derive the dynamical susceptiblity following the 
procedure of ref.ll. Figure 8 shows the two-particle 
spectrum Imx'^iuJ + i6)/ujC at low temperature with 
Ji = — J2. The susceptibilities should satisfy the rela- 
tion Xm ~ Xm ' which have been calculated separately to 
check the accuracy of the Monte Carlo simulation. The 
corresponding results for Xm ^^'^ ^'^^ close enough to 
each other. Hence the validitiy of the Fade approxima- 
tion is confirmed. The sign ImxJ^(ix') < corresponds 
to the negative sign of Rexj|(0) < by the Kramers- 
Kronig relation, and indicates an antiferromagnetic cor- 
relation between the pseudo spins. We find an inelastic 
peak for about Ji < 0.25 (see Fig. 8(a)-(c)). This peak 
corresponds to excitation from the CEF singlet to the 
triplet. On the other hand, no inelastic peak is observed 
at Ji = 0.3 as shown in Fig. 8(d). Instead, the spectrum 
shows only a broad quasi-elastic peak. 

Let us compare the two-particle spectrum with the 
i-matrix, which corresponds to the single-particle exci- 
tation of the local state. As shown in Fig. 7, the single- 
particle spectrum has a gap whose shape depends on the 
strength of the interaction. With r <C 1 as shown in 
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Fig. 7. The imaginary part of the f-matrix at T = 0.001 in the 
region (II) with the condition (a) Ji = — J2, and (b) Ji = 0.3. 
The two Hnes for each (Ji, J2) correspond to up and down spins, 
and they should coincide in the exact results. The parameters 
(Jl, J2) are shown as (a) triangles, and (b) squares in Fig. 1. 



Fig.7(a)(0.1, —0.1) and Fig. 8(a), the single-particle spec- 
trum looks like a step function, and the two-particle spec- 
trum has the inelastic peak corresponding to the gap. 
This means that the ordinary CEF picture is valid, and 
the Kondo effect is insignificant. On the other hand, with 
r - 1 as shown in Fig. 7(a) (0.3, -0.3) and Fig. 8(d), the 
two-particle spectrum has a broad quasi-elastic peak. We 
emphasize that the single-particle spectrum still has a 
clear gap. Namely, the system with r ~ 1 shows both 
features of the Kondo and CEF effects. In this case, the 
simple CEF picture is no longer valid. 

The validity and breakdown of the CEF picture is 
also seen in the components of the static and dynamical 
susceptibilities. In Fig. 8, we have imposed the condition 
Jl = —J2, which leads to It = 0. In the case of r ^ 
1 as shown in Fig. 8(a), combination of Xm according 
to eq. (17) gives small Xt- This is understood because 
the triplet lies higher than the singlet. With stronger 
coupling as shown in Fig. 8(b)-(d), on the other hand, Xt 
is no longler negligible although the interaction contains 
only X^. This increase Xt is due to the Kondo effect. 

Let us turn to the case of Ji < | J2I and Jipc ^ 1, 
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Fig. 8. The two-particle spectrum with the condition Ji = —J2 
at T = 0.001. The parameters (Ji, J2) are shown as triangles in 
Fig. 1. 
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0.001 



Fig. 9. Temperature variation of static susceptibilities with 
(^1,^2) = (0.3, -0.3). 




Fig. 10. The fitted values of A^^, Acef defined by eq.(5), and 
ui/ii, as a function of the coupling constant Ji = — J2. In the CEF 
picture, all these excitation energies should be the same. The 
CEF picture becomes reasonable for Ji < 0.15. The parameters 
are shown as triangles in Fig. 1. 



corresponding to the point with asterisk in Fig. 1. Here 
It is negative and r <C 1. In this case, Xt is small, and the 
inelastic peak can be seen in the two-particle spectrum. 
The behavior of this parameter is similar to the case of 
{Ji, J2) ~ (0.1, —0.1). This similarity is understood since 
a negative It is renormalized to zero at low temperature. 
Therefore, we may say that 7^ is caused by a large 
antiferromagnetic coupling Ji . Owing to this interaction, 
the ground state of the system changes to more stable 
Kondo singlet from the CEF singlet. As in the ordinary 
Kondo system, this singlet ground state in the ST Kondo 
model is strongly coupled with conduction electrons. 

6.3 Static Susceptibility 

The static susceptibility with (Ji, J2) = (0.3,-0.3) 
is shown in Fig. 9. The susceptibility tends to a con- 
stant at low temperature. This is related to the van 
Vleck paramagnetism which arises from virtual transi- 



tions to the triplet CEF level. We can see this param- 
agnetism for other parameters in the region (II) such as 
{Ji,J2) = (0.1,-0.1). Note that the paramagnetic be- 
havior can be seen even in the case where there is no 
inelastic peak in the two-particle spectrum. 

In the singlet-triplet system with the CEF split- 
ting A, the van Vleck susceptibility is given by x = 
2|(s| Jz|tO)p/A, where |s) and |tm) are the singlet state 
and the triplet states (m = -|- 1,0,-1), respectively. To 
study the van Vleck susceptibility in our model more pre- 
cisely, we define the parameter Afj_^, which plays the role 
of effective CEF splitting, as follows: 

-^(s|5^|tO)(tO|^J|s) (22) 



Xm \t= 



In the ST Kondo model, the matrix elements in eq. (22) 
equals to (— 1)'^+''/4. Figure 10 shows the fitted results 
for Af^iy, together with the excitation energy w^j/ corre- 
sponding to the inelastic peak in the two-particle spec- 
trum, and AcEF defined in eq.(5). In the weak coupling 
range with Ji < 0.15, these values are almost the same. 
Hence, in this case, the effective CEF splitting is close to 
AcEF derived by the second-order perturbation theory. 
In other words, the Kondo effect is unimportant. This 
energy has appeared also in the i-niatrix shown in Fig. 
7(a). Namely, the magnitude of the energy gap is nearly 
equal to 2Acef, and the van Vleck paramagnetism ac- 
counts for the static susceptibility. 

As the coupling increases, however, such values as 
A^jy and ujfj,^ deviate from each other. Because of the 
Kondo effect, the CEF picture cannot explain the low- 
temperature susceptibility in this region. 

7. Summary 

In the present paper, we have extended and applied 
the CT-QMC method to the ST Kondo model where the 
conduction electrons interact with two pseudo spins. We 
have derived the impurity t-matrix, spin susceptibilities 
and singlet occupation rate in our model. Since the CT- 
QMC does not use any approximation, the results given 
in this paper are highly reliable with only statistical er- 
rors. 

In the doublet ground-state region (Ji,J2 > 0), 
the behavior of this system is understood as the under- 
screened Kondo effect. In the region J1J2 < for the 
singlet ground-state, the CEF effect competes with the 
Kondo effect. In the case of Tk ^ Acef, the single- 
particle spectrum has a gap structure at low tempera- 
ture, and the two-particle spectrum has an inelastic peak 
corresponding to the gap. In the case of Tk ~ Acef, 
however, the inelastic peak vanishes in the two-particle 
spectrum, or the dynamical susceptibilty. Instead, a 
broad quasi-elastic peak appears due to the Kondo ef- 
fect. The single-particle spectrum retains a gap even in 
this regime. This contrasting behavior between single- 
and two-particle spectra results from the competition be- 
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tween the Kondo effect and the CEF effect. 

Using the algorithm given in this paper, it is possi- 
ble to apply the CT-QMC method to more complicated 
models such as the quadrupolar Kondo model proposed 
by Cox.^^^ As another extension, we can apply the CT- 
QMC to the lattice system using dynamical mean-field 
theory (DMFT), in which the problem is reduced to an 
effective impurity problem. We shall study the ST Kondo 
lattice system in the future work. 
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Appendix: CT-QMC for Ferromagnetic interac- 
tion 

As noted in rcf.ll, the CT-QMC for the Coqbhn- 
Schrieffer model is applicable only to the antiferromag- 
netic coupling (J > 0) because of the negative sign prob- 
lem. In this section, wc show that the algorithm can be 
applied to the ferromagnetic coupling (J < 0) by a slight 
modification. 

The exchange interaction term can be written in the 

following form: 

J^X^^>cl,Ccr = J^Xcra'{cl,Ca- - ad^a') + OiJ 

The constant term may be neglected. Here we have in- 
troduced the parameter a, which gives the equal-time 
Green function as discussed in ref. 9. This parameter is 
set as 

\ (J < 0). 

The two choices give either g{T = +0) for J > or 
g{T = —0) for J < as the equal-time Green function. 
Wc have found that a simulation using this expression 
does not encounter the negative sign configuration for 
N = 2. The absence of the negative sign is understood 
by considering the fc = 1 term in cq.(13). 

We note that the procedure is valid only in the case 
of A'' < 2, where N is the number of the local states. 
Namely, the interaction has the operator X^i and X^-f 
which change the local state. In the N = 2 case, the total 
number of these operators in a configuration q must be 
even. Then the contribution to the weight W{q) has no 
difference between the J > and J < cases. Hence we 
can deal with the ferromagnetic interaction in the case 
of A'' = 2. This fortuitous situation does not occur for 
A" > 3. 
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